Quantum analysis of Rydberg atom cavity detector 
for dark matter axion search 



Akira Kitagawa,* Katsuji Yamamoto^^ 
Department of Nuclear Engineering, Kyoto University, Kyoto 606-8501, Japan 

Seishi Matsuki^ 

Nuclear Science Dtvision, Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan 

Quantum calculations are developed on the dynamical system consisting of the cosmic axions, 
photons and Rydberg atoms which are interacting in the resonant microwave cavity. The time 
evolution is determined for the number of Rydberg atoms in the upper state which are excited 
by absorbing the axion-converted and thermal background photons. The calculations are made, 
, in particular, by taking into account the actual experimental situation such as the motion and 

^\ • uniform distribution of the Rydberg atoms in the incident beam and also the spatial variation of 

0^ \ the electric field in the cavity. Some essential aspects on the axion-photon-atom interaction in 

the resonant cavity are clarified by these detailed calculations. Then, by using these results the 
j^jQ' detection sensitivity of the Rydberg atom cavity detector is estimated properly. This systematic 

, quantum analysis enables us to provide the optimum experimental setup for the dark matter axion 

search with Rydberg atom cavity detector. 
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I. INTRODUCTION 



> 

in 

■<!::;j- ■ Search for the so-called "invisible" axions as non-baryonic dark matter is one of the most challenging issues 
00 in particle physics and cosmology [^,0. The mass range rria ~ 10~^eV — 10~^eV is still open for the cosmic axions 
(The light velocity is set to be c = 1 throughout the present article.) As originally proposed by Sikivie ||], the 
basic idea for the dark matter axion search is to convert axions into microwave photons in a resonant cavity under 
, a strong magnetic field via the Primakoff process. It is, however, extremely difficult to detect the cosmic axions due 
to their unusually weak interactions with ordinary matters. Pioneering tries with amplification-heterodyne-method 
were published already [Q. An advanced experiment by the US group is currently continued, and some results have 
f-^ • been reported, where the KSVZ axion of mass 2.9 x lO^^eV to 3.3 x lO^^eV is excluded at the 90% confidence level 
(]J [ as the dark matter in the halo of our galaxy . 

• In this paper, we describe a quite efficient scheme for the dark matter axion search, where Rydberg atoms are 
' utilized to detect the axion-converted microwave photons [p^|-p^. An experimental apparatus called C ARRACK I 
(Cosmic Axion Research with Rydberg Atoms in resonant Cavities in Kyoto) is now running to search for the dark 
matter axions over a 10 % mass range around lO^^eV. Based on the performance of CARRACK I, a new large- 
scale apparatus CARRACK II has been constructed recently to search for the axions over a wide range of mass . 
Clearly, in order to derive quantitative and rigorous results from the axion search with this type of Rydberg atom 
cavity detectors, it is essential to develop the quantum theoretical formulations and calculations on the interaction 
of cosmic axions with Rydberg atoms via the Primakoff process in the resonant cavity [|l^,^. Here, we present the 
details of these quantum calculations. This theoretical analysis actually provides important guides for the detailed 
design of Rydberg atom cavity detector, which will be examined separately in forthcoming papers. It is also fascinating 
that these quantum theoretical investigations on the axion-photon-atom interaction will be useful in viewpoint of the 
applications of the cavity quantum electrodynamics to fundamental researches. 

This paper is organized as follows. In Sec. ||, we describe the Rydberg atom cavity detector for the dark matter 
axion search. In the quantum theoretical point of view, this Rydberg atom cavity detector can be treated as the 
dynamical system of interacting oscillators with dissipation. In Sec. El|, these quantum oscillators describing the 
photons, axions and atoms are introduced appropriately. Then, in Sec. Hw the interaction Hamiltonians are provided 
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in terms of these quantum oscillators. In Sec. the quantum dynamics of interacting oscillators with dissipation is 
generally formulated in the Liouville picture. A series of master equations is derived to determine the time evolution 
of the quatum averages of the occupation numbers and higher-order correlations of these oscillators. By applying this 
formulation to the axion-photon-atom system in the resonant cavity, we can, in particular, calculate the number of 
the atoms in the upper state which are excited by absorbing the axion-converted and thermal background photons. 
In Sec. VI, we examine some characteristic properties of the axion-photon-atom interaction by sol ving analytically 



the master equation for the simple case with time- independent atom-photon coupling. Then, in Sec. VII, in order to 
make precise estimates on the sensitivity of the Rydberg atom cavity detector, we elaborate the calculations by taking 
into account the motion and uniform distribution of th e Rydberg atoms in the incident beam and also the spatial 
variation of the electric field in the cavity. In Sec. VIII, the dependence of the detection sensitivity on the relevant 



experimental parameters is discussed. Detailed numerical calculations are performed in Sec. IX, and the quantum 
evolution of the axion-photon-atom system in the resonant cavity is determined precisely. The counting rates of signal 
and noise are calculated with the steady solutions of the master equation. The sensitivity of the Rydberg atom cavity 
detector is then estimated from these calculations. Finally, we summarize the present quantum analysis in Sec. 
Appendices are devoted to some supplementary issues. 



II. RYDBERG ATOM CAVITY DETECTOR FOR DARK MATTER AXION SEARCH 



A schematic diagram of the Rydberg atom cavity detector (CARRACK I) is shown in Fig. 0. The axions are first 
converted into photons under a strong magnetic field in the conversion cavity. Then, the photons are transferred to 
the detection cavity via a coupling hole, and they are absorbed there by Rydberg atoms. The detection cavity is set 
to be free from magnetic field to avoid the complexity of the atomic energy levels due to the Zeeman splitting. The 
Rydberg atoms, the transition frequency of which is tuned approximately to the cavity resonant frequency ~ IGHz, 
are prepared initially in a lower state with principal quantum number n ~ 100 by exciting alkaline atoms in the 
ground state with laser excitation just in front of the detection cavity. It should, however, be noted that atoms in 
the upper state with n' (> n) are not generated at this stage. The atoms prepared in this way are excited to the 
upper state by absorbing the microwave photons in the detection cavity, and they are detected quite efficiently with 
the selective field ionization method jl^ just out of the cavity. By cooling the resonant cavity system down to about 
10 mK with a dilution refrigerator in high vacuum, the thermal background photons can be reduced sufficiently to 
obtain a significant signal-to-noise ratio. Hence the Rydberg atom cavity detector, which is free from the amplifier 
noise by itself, is expected to be quite efficient for the dark matter axion search. 

In the quantum theoretical point of view, the Rydberg atom cavity detector can be treated as the dynamical system 
of the interacting oscillators with dissipation which appropriately describe the axions, photons and Rydberg atoms. 
In the following sections we develop detailed quantum theoretical formulations and calculations for this dynamical 
system. Specifically, the analysis is made by taking into account the actual experimental situation such as the motion 
and uniform distribution of the Rydberg atoms in the incident beam and also the spatial variation of the electric 
field in the cavity. This quantum treatment provides proper estimates on the sensitivity of the Rydberg atom cavity 
detector for dark matter axion search. 



III. AXION-PHOTON-ATOM SYSTEM IN THE CAVITY 



We first identify the quantum oscillators which appropriately describe the photons, axions and Rydberg atoms in 
the resonant cavity. This provides the basis for investigating various properties of the axion-photon-atom system. 



A. Resonant mode of photons 

The electric field operator in the cavity V is given by 

E{x, t) = (?icJc/2eo)^/^[a(x)c(0 + a.*{^)c\t)] (3.1) 

for the radiation mode with a resonant frequency where eq is the dielectric constant. The creation and annihilation 
operators satisfy the usual commutation relation 

[c,ct] = l. (3.2) 
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The mode vector field q:(x) is normalized by the condition 



/ |a(x)|2d3x= 1. 
Jv 



(3.3) 



The whole cavity V may be viewed as a combination of two subcavities, the conversion cavity Vi with volume Vi and 
the detection cavity V2 with volume V2, which are coupled together: 

V = Vi©V2. (3.4) 

The axion-photon conversion takes place in Vi under the strong mangetic field, while the Rydberg atoms are excited 
by absorbing the photons in V2. It is then suitable to divide the mode vector as 

a(x) q;i(x) + q;2(x), (3.5) 

where ai(x) = for x G V2 and 0:2 (x) = for x e Vi, respectively. The normalization condition of a{x.) is rewritten 
as 

\ai{x)\^d^x+ [ |a2(x)pd^a; = 1. (3.6) 

Vl JV2 

The actual cavity is designed so that neglecting the small joint region the subcavities Vi and V2 admit the mode 
vectors Q:i(x) and 0:2 (x) (up to the normalization and complex phase), respectively, whose frequencies are tuned to 
be almost equal to some common value In this situation, as confirmed by numerical calculations and experimental 
observations, two nearby eigenmodes with the frequencies luc, ^'c — '^c obtained for the whole cavity V. Then, the 
mode vector a(x) is constructed approximately of ai(x) ~ cki(x) and a2(x) ~ ck2(^) with significant magnitudes 
in both Vl and V2- The conversion of the cosmic axions takes place predominantly to the radiation mode which is 
resonant with the axions satisfying the condition \ujc~fna/h\ < 7^ (axion width) small fraction of 7 (cavity damping 
rate), as will be discussed later. The cavity can be designed so as to give a sufficient separation of jwc — o.'cl > several 7 
for the nearby modes with strong coupling between Vi and V2. Therefore, in the search for the signal from the cosmic 
axions, the one resonant mode can be extracted solely for the electric field in a good approximation, as given in Eq. 



(3.1), whose frequency ujc is supposed to be close enough to the axion frequency uja — nia/h. 



B. Coherent mode of cosmic axions 



The axion field operator is expanded as usual in terms of the continuous modes: 



(x,0 = h^/^ 



0+(x, t)+<f)-{x,t), 



ak(t)e' 



(3.7) 



where cf)'^ and (f> represent the positive and negative frequency parts, respectively. The cosmic axions form a coherent 
state with momentum distribution ?7a(k). 



k>=iV-/.exp(/^^,.(k)4) 10), 



(3.8) 



where iV,, is the normalization constant to ensure the condition (?7a|'7a) = 1- The velocity dispersion of the galactic 
axions is expected to be very small as Pa ~ 10~^ so that f7a(k) takes significant values only in a small region TZa where 

\k\<(3afna/n m. 



The axion field operator (3.7) may be divided into the low- momentum part (j)fi^ with k G TZa and the residual part 



(x, t) = (x, t) + (/)i.cs(x, t). 



(3.9) 



Since the coherent region TZa is chosen so as to give (?7a|</>|?7a) ~ (Val't'TialVa) , the residual part ^res does not provide 
significant contributions in the following calculations, i.e., (?/a|0rcs|?7a) ~ 0. It is further noticed that the spatial 
variation of 07j_^(x, i) is negligible in the cavity region, i.e., e*'' '' ~ 1 for |k| < (3ama/fi and x e Vi. This is because 
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the de Broglie wavelength of the axions, Xa — h/{f3a'ma) ~ lOOm typically for iria ^ 10~^eV and f3a ^ 10~^, is much 
longer than the microwave cavity scale 0.1m. In this situation to ensure approximately 07j^ (x, t) « (p-ji^ (0, t) in the 
cavity, the coherent axion mode can be identified as 



a{t) = (?iI],)-i/2^+JO,i) 

r dPk 



^a'" I T^^^Mty (3.10) 



Here, the normalization factor is given by 



„ I (Pk 1 / B„m„ \ , , 



so that the coherent mode operator satisfies the canonical commutaion relation, 

[a,at] = l. (3.12) 

The normalization factor is expressed in terms of the mean velocity (3a of galactic axions with Uk — rua/h and 
f3a ^ 1. It should, however, be realized that this parameter f3a is introduced just for the normalization of the coherent 
axion mode. The actual mean velocity of axions is rather determined by the axion spectrum, which may slightly be 
different from the parameter Pa given in Eq. ( |3.11 ) . The normalization factor is canceled out anyway in calculating 
the axion signal rate, as explicitly seen later. 

The coherent cosmic axions can be described effectively in terms of a single mode oscillator, as shown in the above. 
Then, the energy spread of the galactic axions can be taken into account as the damping rate of the coherent axion 
mode, 

7a - Plrria/h. (3.13) 

The expectation value of the number operator a^a of the coherent axion mode is determined by the energy density pa of 
the cosmic axions as follows. The number density operator for the axion field is given by na(x, t) — {2h)~-^ {i(f>^ docf)'^ — 
ido4>~4'^), and its expectation value is calculated for x S Vi and /3a <C 1 as 

{ria\na{-x,0)\'qa) ^ h^'^ma{r]a\(t>^^{O,0)(l)^^{O,0)\r]a) 

^ Pa/ma. (3.14) 



Then, by considering Eq.( 3.10 ) to relate the relevant part of axion field operator (piz^ to the coherent mode operators 
a and , we can calculate the occupation number of the coherent axion mode as 

na - {va\aHO)amrja) ( ( —) ■ (3-15) 



PamaJ \m. 

This expression implies that the coherent axion mode is normalized suitably in a box with a volume of 
(de Broglie wavelength)"^. 



C. Oscillator for Rydberg atoms 



The Rydberg atom is treated well as a two-level system in the resonant cavity. The relevant lower and upper 
atomic states are represented by \n) and respectively. They are connected by the electric dipole transition with 
frequency LUt — (En' — En)/h, which is actually fine-tuned to be almost equal to the cavity frequency tOc by utilizing 
the small Stark shift. This two- level atomic system is described in terms of spin 1/2 like operators 

D+^\n'){n\, D-^\n){n'\, 

D'^^i\n'){n'\~\n){n\). (3.16) 

These operators satisfy the usual SU(2) Lie algebra, 

[D+,D-] = 2D^, [D^ D^] = ±D^. (3.17) 
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The free Hamiltonian for the two-level system is then given by 



(3.18) 



The Rydberg atoms are initially prepared in the lower state in the present detectoin scheme for cosmic axions. It 
is also expected that the probability of the excitation of an atom to the upper state by absorbing a photon is small 
enough. This is because the average number of photons in the cavity is much smaller than one at sufficiently low 
temperatures (~ lOmK). Therefore, the atoms almost remain in the lower state, providing a good approximation 



(3.19) 



in the first commutation relation of Eq.(3.17). Then, the two-level atomic system operators may be substituted by 
those of an oscillator as lITsIl 



which satisfy the commutation relation 



(3.20) 



(3.21) 



The second commutation relation of Eq. ( 3.17 ) is also reproduced with 
and accordingly the free atomic Hamiltonian is expressed effectively as 



Tiujh 



bW up to the higher order terms. 



(3.22) 



where the irrelevant constant Tiujh/'i- is discarded. This treatment of the Rydberg atom in terms of the quantum 
oscillator is valid if the probability in the upper state (or the lower state) is very small, i.e., {b'^b) ^ 1, as is the case 
in the present scheme. 



D. Characteristics of the photons, axions and atoms 



We have seen so far that the resonant photons, coherent cosmic axions and Rydberg atoms in the resonant cavity 
are described in terms of the appropriate quantum oscillators with dissipation due to the couplings to the relevant 
reservoirs. The characteristic properties of these quantum oscillators are summarized as follows. 

The thermal photon number fic of the resonant mode is determined by the cavity temperature Tc as 



The damping rate 7c of photons is described in terms of the quality factor Q of the cavity as 



The coherent axion mode is normalized in a box of the de Broglie wavelength, as seen in Eq. (B.15). The axion 
number is then estimated as 

where the enegy density of the cosmic axions Pa is taken to be equal to that of the galactic dark halo phaio — 
O.SGeVcm^''. The dissipation of the coherent axions is characterized by their energy spread as 

7a ^ Plmjh = 0.027 (—J^^S__\ , (3.26) 



^10-6 X 2 X 10^^ 

where the resonant condition ~ TilJc is considered. All the atoms are prepared in the lower state so that 

m = 0. (3.27) 
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The atomic dissipation is determined by its lifetime Tb as 7fc = . The Hfetime of the Rydberg state is typically 
Tb ~ lO^^s for n ~ 100 in the vacuum [17|. Since the transitions to the off-resonant states are highly suppressed in 
the resonant cavity, the atomic lifetime may even be longer. Hence, the damping rate of Rydberg atoms is estimated 
to be at most 

76 = 6.6 X lO^^eVn^M — , (3.28) 
which is actually much smaller than the photon damping rate 7. 



IV. INTERACTIONS IN THE RESONANT CAVITY 



In this section, we provide the interaction Hamiltonians in terms of the quantum oscillators. They determine the 
quantum evolution of the axion-photon-atom system in the resonant cavity, which will be investigated in detail in the 
following sections. 



A. Axion-photon interaction 



The axion-photon interaction under a strong static magnetic field with flux density Bq is described by the La- 
grangian density 



Ca = ?i^^^eo5a77'/'-E • Bo 



(4.1) 



where Ti'^^ and eo (dielectric constant) are explicitly factored out so that the Lagrangian density has the right 
dimension Ca ~ ?is^^m^'^ ^ eVm^'^ with the axion-photon-photon coupling constant ^077 ~ eV~^, as given below. 
The axion-photon-photon coupling constant is calculated P,p| as 



_ a nia {1 + Z) 

?a77-Ca77 2^^^^^ ^ , 



where Z = mu/nid, and 



with 



E 2{A + Z) 
C ^ 3(1 + Z) 



E ~ TrQpqQpj^, CSab = TrQpqAaAfc 



(4.2) 



(4.3) 



(4.4) 



The parameter Ca-yy represents the variation of the axion-photon-photon coupling depending on the respective Peccei- 
Quinn models such as the so-called KSVZ g and DFSZ [|. 

The original Lagrangian density for the axion-photon-photon coupling in Eq. (4.1) provides the effective interaction 
Hamiltonian between the coherent axion mode a and the resonant radiation mode c. 



Hac = hK{a^c + ac^). 



(4.5) 



The axion-photon conversion in the cavity Vi is well described with this interaction Hamiltonian. The coupling 
constant n is determined for ujc — raa/fi by considering the relations |(?7a|0^|'7a)| — ^(Pa/'Ti^)"'^/^ from Eq. (3.14) and 
\{'na\a\'na)\ = nj from Eq. ( |3.15| ) in calculating {ria\ J^^{-Ca)d^x\ria) 



_ tl/2 1/2 n 

K — n 9ajj^Q ^cff 

= 4 X 10"26gY^-i 



f3ama \ V\ 
2TTh / 2 



1/2 



9a 



77 



10-3 X lO-^cV 



1.4 X 10-i5GeV"^ 
3/2. y/2 

V5000cm3 J 



4T 



(4.6) 
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where 



Bcff — CiGBq 



(4.7) 



and Bq is the maximal density of the external magnetic flux. The axion-photon-photon coupling constant ga-yy is 
taken here to be the value expected from the DFSZ axion model Q at iria = 10~^eV. 

This effect ive axion-photon coupling k apparently involves the (3a dependence comin g from the normalization given 
in Eq. (3.11). This Pa dependence is, however, canceled with that of fia given in Eq. (3.15) in calculating the signal 



rate which is actually proportional to (K/7) tIq. The form factor for the magnetic field is given by 



/ d3a;c,i(x) • [Bo(x)/So] 



with 



Ci 



/ d^x\a.i{yL) 
Jvi 



-,1/2 



(4.8) 



(4.9) 



This additional factor Ci (< 1 as seen from Eq.( |3.6| )) represents the effective reduction of the axion-photon conversion 
which is due to the fact that the magnetic field is applied only in the conversion cavity Vi . We may obtain the effective 
magnetic field strength i?eff — 4T, as taken in Eq.(^|^), by using typically a magnet of Eq ~ 7T and the cavity system 

with G = \/0.7 of TMqio mode and Ci — 0.7 for the conversion cavity. 

It should here be remarked that the possible interaction term hK,{ac + a^c^) is not included in Eq. (4.5). This is 
justified as follows. The characteristic time scale for the evolution of the system considered in the present scheme 
is placed by the lifetime of photons — 7"^ in the cavity. Hence, although the interaction term hK{ac + a^c^) 
representing the processes with energy change of huJc + rua is also obtained from the original axion-photon-photon 
coupling Ca, its effects are suppressed sufficiently for Q 3> 1 by the energy conservation in this time scale r^. 



B. Atom-photon interaction 



We next consider the interaction between the Rydberg atoms and the resonant photons. The Rydberg atoms are 
utilized for counting photons in the cavity. They are excited by absorbing the photons through the electric dipole 
transition with frequency Ub- The emission and absorption of photon by the two-level atomic system is described by 
the interaction Hamiltonian 



Hn = nn[D+c + D-c^] 

~ hn[b^c + bc'']. 



(4.10) 



Here, we may assume for simplicity, though not essential, that the mode vector field is described in terms of a uniform 
complex polarization vector €2 (|e2| = 1) and a real profile function /2(x) as 



q;2(x) = esKj ^^V2(x) 



(4.11) 



with V2 = |Q!2(x)|j7iax- The profile function is normalized for convenience by the condition |/2(x)|inax = 1 at the 
antinodal positions. Then, the intrinsic atom-photon coupling constant Q is evaluated at the antinodal position in 
terms of the electric dipole transition matrix element d = \e2 ■ d\ projected to the direction of polarization: 



1/2 



where V2 ~ V2 is expected in accordance with the electric field normalization condition (3.6) 
coupling at an arbitrary atomic position in the cavity is also given with the profile function by 

r!(x) = r!/2(x). 



(4.12) 



The atom-photon 



(4.13) 



In the actual experimental system, the Rydberg atoms are injected as a uniform beam. Then, certain number N 
of Rydberg atoms are constantly present in the cavity. Such an ensemble of identical atoms behaves as a collective 
system in the interaction with the resonant radiation mode. The novelty of Rydberg atom physics is to provide 
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situations where this collective behavior appears for a relatively small number of atoms, typically N 10^ — 10^ [ p^ . 
Suppose for simplicity that the N atoms are at the antinodal position (0(x) = Q) in the detection cavity. Then, the 
effective coupling between the radiation mode and the A'' atoms is actually given by 



N 



Hbc = J2 Hnih,c) = nnNib^c + bc^), 



(4.14) 



where the collective atomic mode operator is defined by 



1 ^ 



(4.15) 



satisfying the commutation relation [6,6^] — 1. (Hereafter we may use without confusion the operators b and b^ for 
the collective atomic mode rather than the single atomic mode.) The collective atom-photon coupling = ^^Vn 
[n8| is estimated typically as 



l^AT = 1 X 10"^°eV;i 



5 X lO^s-i 



N 

103 



1/2 



(4.16) 



Hence, if the number of Rydberg ato ms is as large as '--^ 10'^, this collective coupling can be comparable to the 
cavity damping rate 7 as given in Eq. (3.24). When the atomic motion and distribution in the ca vity are taken into 
account, the atom-photon coupling should be modified suitably. This point will be treated in Sec. VH. 



V. QUANTUM EVOLUTION OF THE SYSTEM 



We here formulate the quantum dynamics of interacting oscillators with dissipation. A systematic procedure 
is developed in the Liouville picture to calculate the quantum averages of the particle numbers and higher-order 
correlations. In practice, a series of master equations is derived for such quantities from the Fokker-Planck equation 
based on the Liouville picture. (The quantum evolution of the system can also be described in the Langevin picture, 
which is considered in Appendix ^.) By applying this formulation, the quatum evolution of the axion-photon-atom 
system in the resonant cavity is determined with the effective interaction Hamiltonians presented in the proceeding 
section. Then, we can calculate, in particular, the number of the atoms in the upper state which are excited by 
absorbing the axion-converted and thermal photons. 



A. Evolution in Liouville picture 



The reduced density matrix of the damped oscillators qi (such as a, 5, c) interacting each other obeys the Liouville 
equation |19|, 



dp 1 



[H,p]+Ap. 



dt ih 

The Liouvillian relaxations are represented by the operator Ap, which may explicitly be given by 



(5.1) 



Ap = ^ y ['2q^pql - qlq'i.p - pqlq^ 



E 



(5.2) 



where 7^ and rii are the damping rates and equilibrium occupation numbers, respectively. The total Hamiltonian is 
given by 



H = ^ hujiqlq, + ^ hn,j{t)qlqj, 



(5.3) 



8 



where ^ij{t) = represent the interaction terms, which may be time-dependent, as is the case for the interaction 

between photons and moving atoms in the cavity (see Eq. ( 5.10|) ). 

The solution of the Liouville equation is generally obtained in terms of the creation and annihilation operators as 

p{q^,q\,t). (5.4) 

The quantum average of a physical quantity represented by a relevant operator O is evaluated with this density matrix 
by the formula 

d{t) = {O) = TT[0{q,,q\,t)p{q,,qlt)]. (5.5) 

In practical calculations, it may be useful to take the coherent state basis Then, the density operator p is 

represented by the "classical" time-dependent distribution function P{ai, a*, i), and the quantum average is calculated 
by 

0{t) = [W<faMa,,a*,t)P{a„a*,t). (5.6) 

i 

Here the subscript "n" means the the expression of the operator O when it is written in the normal ordered form. 
For example, On{a, a*) = a*a + 1 for 0{q, gt) = ggt 

The Liouville equation is expressed in the coherent state basis as the Fokker-Planck equation to determine the 
distribution function P: 

zW.,(t)7^(a,P) - tH:^{t)j^{a;P) 



dt '■"'da,:-' ' ''"da 

d^P 
' daida 



Here, the effective Hamiltonian and diffusion term for the damped oscillators are represented by the matrices 

n^J{t) = 1^, - ^7,^ % + n,,{t)(i - 5,j), (5.8) 

'Dij = -yifiiSij. (5.9) 
They are given explicitly for the axion-photon-atom system under consideration by 

n{t) = I n{t) We - 57c n I (5.10) 








with uja = raa/fi, and 



Ibnb 

V=\ 7enc I . (5.11) 

7ana 

The collective atom-photon coupling becomes time-dependent through the atomic motion along the x axis with 
velocity v as 

n{t) = VLN.f{vt). (5.12) 

Here, it is assumed for simplicity that the N atoms are injected together as a pulsed beam. This time-dependence is 
determined by the atomic velocity and the electric field profile f{x) = f2{x, yo, zo) of the detection cavity along the 
atomic beam which is injected in the x direction from the point (0, yo, zq) in the y-z plane. 
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B. Master equations 

We may be interested in the self-adjoint multiple moments of the oscillators, 

= y Y[d'^o:,a*^ ■ ■ ■ a*^aj, ■ ■ ■ aj^P{a^,a* ,t), (5.13) 

i 

where 2p denotes the number of involved operators. A series of master equations is then obtained for these moments 
by considering the Fokker-Planck equation for the distribution function: 

^^^^ = -i7^(2p)(;)_^(2p) ^ p(2p-2)_^(2p-2)^ (5^^4) 

dt 

Here the linear operators Ti^^^^ and I?(2p-2) a,cting on the moments are defined by 

k 

-En.;W-^?^V........> (5-15) 

k 

(p(2p-2)_^(2p-2))^^ EE Vl?,,,^^'^'^ (5-16) 

k,l 

with abbreviation of the indices I = ii ■ ■ ■ ip, J = ji ■ ■ ■ jp, = ii ■ ■ ■ ik-iik+i ■ ■ ■ ip and J^' = ji • • • ji^iji+i ■ ■ ■ jp- 

We examine, in particular, the second-order moment A/ij (t) = A/'j^ (t). The number of the atoms which are excited 
by absorbing the photons is given by the diagonal component of Af{t) as 

nt,{t) = (b%) = Afbbit). (5.17) 

The master equation for Af{t) reads 

^ = -tJ\fH^{t)+tH*{t)J\f + V. (5.18) 
at 

This linear equation with an inhomogeneous term may be solved formally as follows. We first introduce a new matrix 
Af'{t) instead of J\f{t) by 

M'{t)^U-^*{t)N{t)U-^^{t). (5.19) 
The linear transformation tl{t) representing the time evolution due to 7i(t) is given by 



Uif) = P 



— zcxp / H{T)dT 

Jo 



(5.20) 



satisfying the equation dU/dt = —iH(t)U with U{0) = 1, where "P" denotes the chronological product. Then, the 
master equation ( 5.18| ) is reduced as 

^=U-^*{t)VU-^^{t). (5.21) 

By considering the initial condition Af'{Q) ~ A/'(0) with U{Q) — \ and Eq. ( |5.9| ) for V, we can solve the above equation 
as 

Ml^{t)^B%{t)nu+N.M. (5.22) 

where 

S,^-(t)= [\uUg'*{T)U^l^{T)dT. (5.23) 
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Then, we obtain the solution as 



where 



n%[t)^[W[t)B\t)U'^{t)]^ 



(5.24) 



(5.25) 



The same result is also obtained in the Langevin picture (see Appendix This solution approaches to certain 
asymptotic value after a long enough time as 



M,(i) « T^Muk {t » max[7-i]). 



(5.26) 



In practical calculations, the master equation (5.18) will be solved numerically with the time-dependent effective 
Hamiltonian 7i(t) as given in Eq. ( 5.1(]| ). 



VI. ASPECTS OF AXION-PHOTON-ATOM INTERACTION 



We can see some characteristic properties of the axion-photon-atom interaction in the resonant cavity by examining 
t he si mple case with the constant atom-photon coupling ^(t) ^ for the effective Hamiltonian Hit) = H in Eq. 
( 5.10 ). In this case, as derived in the Appendix the analytic solution is obtained for the particle numbers (with 
the condition = 0) as 



(6.1) 



where 



1 - 



7j 

Am7 



7j 



A„ 



(6.2) 



with 



= lim {s + iXk)isl+in)^.^, 



Amn — ^"i-i^rn ^ ^n) ■ 



(6.3) 
(6.4) 



Here, the atomic damping rate 7^ may be neglected for simplicity, since it is sufficiently smaller than 7^ and Qj^ 
(76 ~ O.OOI7 with Tf, ^ lO^'^s for ujc ^ 10~^eV and Q ~ 10'*). The condition ujb = lUc rnay also be taken for 
definiteness, since the atomic transition frequency should be tuned almost equal to the cavity frequency. Then, the 
eigenvalues of the Hamiltonian Ti. are given by 



Ai = ujc 

A3 = (^a 



4 ' 



(6.5) 
(6.6) 
(6.7) 



where 



''^'^^ { i^Nh > 1/4) ■ 



(6.8) 



If the number of atoms A'' (or the atomic beam intensity /Ryd) is not so large giving ^n/i < 1/4, the damping rate 
of the eigenmode of Ai is smaller than 7/2, and that of A2 lies between 7/2 and 7. On the other hand, in the strong 
coupling region of SI at /7 > 1/4 for the collective atom-photon interaction the eignemodes of Ai and A2 form a doublet 
around the frequency uJc with the same damping rate 7/2 (Rabi splitting). 

Among the damping rates Re[Amn] of the respective terms in the factors Tij{t) representing the contributions of 
thermal photons and axions, Re[A33] — 7^ (~ O.OI7 for (ia ~ 10^^ and Q ^ 10^, typically) is the smallest one. The 
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rate Re[Aii] ~ /"iY for ^n/i ~ 0.1 may also be comparable to the smallest rate Re[A33]. The atomic transit 

time through the cavity is, on the other hand, given by 

itr = L/v (6.9) 

with the detection cavity length L and the atomic velocity v. This transit time provides the effective cut-off for 
the axion-photon-atom interaction in the cavity. (Here we assume for simplicity that the atoms have the uniform 
velocity.) It is typically itr — 400r^ with L = 0.2m and v = 350ms~^ for rua = 10~^eV and Q = 2x 10'' in the case of 
the detection apparatus such as CARRACK I. Hence, the transit time can be regareded to be long enough compared 
to (Re[A„„])-\ i.e., 

ttr > several 7""', (6.10) 
so that the respective particle numbers will almost reach the asymptotic values as 

r.Attr) « r,,(oo) = ^29^9"^,^^- (6-11) 



By using Eqs. ( |6.5[ ) - (6.7) and the explicit matrix form (A17) for (si + ili.) ^ given in Appendix we can 



calculate the coefficients g^j in Eq. (|6.3D . Then, we can show the relations 

ndoo) = r-cc(oo) = 1. (6.12) 

This implies that if the axion-photon interaction is turned off, the numbers of the photons and excited atoms reach 
the same asymptotic value fic of the thermal photon number at Tc p^ : 

ni,[c b] ^ nc[c ^ c] « fic. (6.13) 
The number of axion-converted photons is, on the other hand, given approximately by 

nc[a ^ c\k rca{oo)na- (6-14) 
The number of excited atoms due to the axion-converted photons is also given by 

nb[a ^ c ^ b] fii rba{oo)na. (6.15) 
These contributions from the axions are essentially determined by the factors 

C C « TT T^T TT' ^^-16) 

where m ^ k,l and k ^ I. The detuning of axion mass from the cavity frequency is given by 

Aa;a = Wa-Wc- (6-17) 

Then, the above factors are enhanced if the axion detuning Awq lies in the range where the resonant conditions 
A3 ~ Al and/or A3 ~ A2 are satisfied. 

In the leading order of the axion-photon coupling k with Eq. ( 3.16| ), the factors for the axion contrib ution s are 
given with rcQ(oo), r;,a(oo) (x (k/j)-^. Then, by noting the relation {K/j)'^ha oc (/5a/wa)Vi from Eqs. ( 3.15| ) and 
(4.6) it is found that the axion-converted photons nc[a c] and the number of atoms nb[a — c ^ &] excited by such 
photons are both proportional to the number of axions contained in the conversion cavity. It is here relevant to define 
the form factors for the axion contributions with respect to the axion detuning (^tr > several 7""^ 3> 7~^) by 



?-m(^tr) ria{oo) 



These form factors ^^^(Awa) for a ^ c — > 6 (solid lines) and aca{^<jJa) for a — > c (dotted lines) are plotted together 
in Fig. ^ for some typical values of the atom-photon couphng, Q.n h = 0.1,0.5, 1.0. 

The behavior of (Jha{^^a) is, in particular, understood by noting its dominant contribution which is in fact given 
by the term of \gZ^^\^ (7a < 7) as 
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aUAu;^) 0. ^^^^^^^ _ ^^^^ ^^^2^^^ _ ^ ^^^^^ . (6.19) 



The peak of this form factor appears at 



(f7jv < H) 



Aa..(peak) = j ±^^^^3^ (^^^ > ^1) ' ^^'^O) 

where = j /y/8 + 0(7q). (Although ^^^(O) is apparently divergent for Qn — 5(707)^^^ in the approximate formula 
of Eq. ( |6.19[ ), the other contributions also become significant around AtOa = so as to give a finite value of crba{0)-) 
It is found that for flj^ ?J 7/4 the peak value of the axion signal decreases due to the Rabi splitting, which approaches 
to crf,a(Aa;a(peak)) ~ 1 for (flisi/j)'^ ^ 1. On the other hand, for the appropriate atom-photon coupling such as 

- (7a7)'^' (6.21) 

the signal is enhanced significantly, as observed in Fig. ^ (iljv = O.I7 and 7a = O.O27), by virtue of the narrow width 
(long coherence) of the galactic axions as 

abaiAuJa) - (7a/7)"\ (6-22) 

when the axion detuning becomes small enough to satisfy the condition 

|Ac^,|<7,. (6.23) 

Here, the energy uncertainty ~ h/tti due to the atomic motion is assumed to be smaller than the axion width h'ya- 

The form factor aca{AuJa) determines the equilibrium number rca{oo)na — 4,acaiALUa){K/j)'^na of the axion- 
converted photons in the cavity. As seen in Fig. ^ its peak value is CTco — 1 almost independently of the atom-photon 
coupling Qn. It is here, in particular, interesting to observe a narrow dip in aca{AuJa) around AtOa = for /l — 0.1. 
This indicates that the axion-converted photons are efficiently absorbed by the atoms for H ^ {la/lY^"^- For 
larger ri^r, two separate peaks appear in (Jca{ALJa) due to the Rabi splitting. 

Some characteristic features concerning the axion-photon-atom interaction in the resonant cavity have been dis- 
cussed so far by making the analytic calculations for the simple case with the constant atom-photon coupling ^Im- 



They will indeed be confirmed in Sec. [X by performing detailed numerical calculations for the realistic situation with 



the continuous atomic beam passing through the spatially varying electric field. 



VII. DETECTION SENSITIVITY WITH CONTINUOUS ATOMIC BEAM 



In order to make precise estimates for the sensitivity of the Rydberg atom cavity detector, we have to take into 
accout (i) the motion and (ii) the almost uniform distribution of the atoms in the incident beam as well as (iii) the 
spatial variation of the electric field in the cavity. We here elaborate the calculations presented in the preceeding 
sections by treating these points appropriately. 

The electric field felt by the atoms varies with time through the atomic motion. Accordingly, the effect of atomic 
motion can be incorporated by introducing the relevant time-dependenc e for the atom-photon coupling, which is 
determined by the profile of the electric field in the cavity, as given in Eq. ( 5.12 ). On the other hand, in order to treat 
the spatial distribution of the atoms, we divide the atoms in the cavity into K bunches with a fixed beam intensity 



-^Ryd = N/tf[ 



(7.1) 



Here, N is the total number of Rydberg atoms in the cavity. Then, the continuous atomic beam will actually be 
realized for if ^ 1. 

The collective atomic mode of each bunch is denoted by 6^ (i = 1, 2, . . . , K), and the effective Hamiltonian is given 
by a {K -I- 2) x (if + 2) matrix 



n{t) = 



f 




^x(t) 


{^h - 


f7b)l 








^K{t) 


f71 ... 




LOc - 57c K 


... 





K bJa- 57a , 



(7.2) 



13 



Now suppose that the i-th atomic bunch locates around 

Xi = {i — l)5x 

for a time interval 

M5t <t<{M + l)6t, 

where M = 0, 1, 2, . . ., and 

5x = L/K, 5t = ttr/K. 



(7.3) 
(7.4) 
(7.5) 



Then, the collective atom-photon coupling for the i-th bunch containing N/K atoms is given with the electric field 
profile by 



where 



a(t) = {Um/^K) f{x,+vt), 
t = t-M5t {0<t< 6t). 



(7.6) 
(7.7) 



After each time interval of St, the A'-th atomic bunch leaves the cavity, and a new one comes in. Then, the collective 
atomic modes of the respective bunches should be replaced as 



(7.8) 



which is represented by a {K + 2) x [K + 2) matrix 



V- = 



/ 








• 


• 











\ 




1 





• 


• 
















1 


• 


• 





















• 


• 





















• 


• 1 





















• 


• 





1 







V 








• 


• 








1 


/ 



(7.9) 



In accordance with this replacement of the collective atomic modes at each period of St, the time evolution of the 
particle number matrix J\fij (t) should be determined by dividing it into the corresponding parts in time as 



J\f{t, M) = JV{t) {M6t < t < {M + l)6t). 

These parts are connected together at t — MSt with a suitable matching condition 

Sf{0, M + 1) = V-AfiSt, M)Vi. 

This condition is written down explicitly as 

Mj(0, M + 1) = JV,-i j-i{St, M) [2 < ^, J < K], 

Mj(0,M + l) =0 [l<]<K + 2], 

Mi(0,M+l) = 0_ [l<i<K + 2], 

Afijlo, M + 1) = M-i j{St, M) [2<i<K,j = a, c], 

Mj(0,M + 1) =M j-i{St,M) [i = a,c,2<j<K], 

Af^-i (0, M + 1) = {St, M) [i, i = a, c] . 



(7.10) 
(7.11) 



(7.12) 



Here, the second and third lines imply that the incoming atomic mode (6i) does not have any correlation initially 
with the other quantum modes (&2, • ■ • , ^x, c, a ) alre ady interacting in the cavity. 

In this way, by solving the master equation (5.18) for many time intervals, we obtain the steady solution 



A/'(i) wA/'(f,A/> 1). 



(7.13) 
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It in fact appears independently of the choice of the initial value Af{0), which is due to the dissipation of the quantum 
modes for the axions, photons and atoms. The number of excited atoms contained in each atomic bunch is then given 
by 

rib, (t) = (t) = nl (t) + nl (t). (7.14) 
Here, the contributions of the axions and thermal photons are proporional to the respective particle numbers as 

= n,a{t)na, nlit) = n^c{t)nc- (7.15) 
The distribution of excited atoms in the cavity is determined with this steady solution as 

p^{x) = ^bAt = x/v) _ 1)^^ < 2. < ^^3.)^ (7 16) 

ox 

The K-th atomic bunch exits the cavity at t = St, and the excited atoms contained there are detected. Accordingly, 
the counting rates for the contributions of the axions and thermal photons are calculated, respectively, for large 
enough K and M by 

Rs = (7.17) 



< (St) 

Rn = (7-18) 

By using these counting rates, the measurement time required to search for the axion signal at the confidence level of 
ma is estimated as 

At= ^ — ^ ' ' . (7.19) 

In the search for the axions with unknown mass, the cavity frequency (wf, = lOc for definiteness) should be changed 
with appropriate scanning step Acjc- The total scanning time over a 10% frequency range is then given by 

itot = -^At- (7.20) 

The frequency step is determined by considering the resonant condition for the absorption of the axion-converted 
photons by the Rydberg atoms. Specifically, as seen in Eqs. ( |6.2l[ ), (|6.22 ) and ( 6.23| ), the signal rate is enhanced 
significantly for the axion detuning |Aa;a| 7a with H ^ (To/t)"'"^'^- Hence, the scanning frequency step should 
be taken as 

AWc ~ la. (7.21) 

Then, if the axion really exists in a mass range searched with this frequency step, the resonant condition for the axion 
signal can be satisfied at a certain scanning step. 



VIII. DEPENDENCE ON THE RELEVANT PARAMETERS 



We here consider how the counting rates of the axion signal and thermal photon noise depend on the relevant 
experimental parameters, before presenting detailed numerical calculations in the next section. Although the notation 
for the case of if = 1 with constant atom-photon coupling VIm is used for simplicity, the essential features argued 
below are indeed valid for the realistic case with continuous atomic beam. 

The energy scales involved in the effective Hamiltonian arc as follows: 

detunings : AcJq = uja — uJc, Awt = Wfc — Wc, 

couplings : k, 

dampings : 7^, 76, 7^=7- 
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The atomic transit time is also important to determine the time evolution of the system, which is controlled by the 
atomic velocity v with a fixed cavity length L as 



V - L/tt 



(8.1) 



The effects of the dissipation of axions and atoms are rather small compared to that of the photons in the actual 
situation with 7a, 7h <C 7. The atomic transit time is in fact much longer than the other characteristic time scales 
except for e.g., itr — 400t^ for v — 350ms^^ with L = 0.2m, Q = 2 x 10"^ and ma = lO^'^eV. In this situation, 
the quantum averaged particle numbers approach to the asymptotic values with the factors (itr) 



, (00), as given 

in Eq. (6.11). Then, the counting rates of the signal and noise, which are calculated by Rs = ?'6a(itr)^a/itr and 
Rn ~ fbc{ttr)nc/tti, are found to be roughly proportional to the atomic velocity v. 

In the weak atom-photon couplin g regi on of Qn/j ^ 0.1, the counting rates Rs and i?„ are roughly proportional 



to (X N (X /Ryd, as seen in Eq. (6.19) for the axion signal. On the other hand, when the atom-photon coupling 



becomes significant with sufficient atomic beam inten sity / Ryd, t he no ise rate saturates to certain value i?„ ~ n-c/^tr, 
and the signal rate Rg is maximized, as seen in Eqs. ( |6.21 ) and ( 6.22 ), for the atom-photon coupling 



f^jv/7 ~ (7a/7)'^' ~ 0.1 



Pa 

10-- 



Q_ 

104 



1/2 



■2) 



It is possible to attain this optimal coupling in the actual detection apparatus by preparing the appropriate number 
of Rydberg atoms with a suitable laser system. The beam intensity of Rydberg atoms ^Ryd = ^ /ttr providing 
this optimal value ( |8.2| ) for the collective atom- photon coupling Vtjsi — fl^/N is chosen depending on the relevant 
parameters as 



/Ryd ~ v(3lm,lQ-'n- 



(8.3) 



It is typically estimated as I^yd ^ 7 x lO^s^^ with N ~ 400 atoms for = 10"^eV, Pa = 10"^, Q = 2 x lO", 
= 5 X lO'^s^^, V — 350ms~^ and L = 0.2m {tti ~ 400T-y). On the other hand, if the atom-photon coupling becomes 
too strong, the signal rate decreases due to the appearance of Rabi splitting, as seen in Fig. ||. 

We have observed that the atom-photon interaction can not be treated perturbatively for VLn/"! ^ 0.1 The 
conversion between atoms and photons appears to be reversible in this case. On the other hand, since the axion-photon 
coupling is extremely small, the conversion of axions to photons can well be treated perturbatively as an irreversible 
process. The signal rate Rs is then found to be proportional to [K/^Y'fia in a very good approximation. Then, by 
considering the relations (3.15), (3.24) and (4.6) with lOc — ma/Ti, the dependence of the signal rate on the relevant 
parameters is specified as 



Rs 



Vci^-,BtsQVl{Palma)(Jba- 



(8.4) 



Here, the form factor aba = o'f,a(±AtJc/2) for Awa — ±Att'c/2 (Awf, — 0) is also indicated. This form factor ac tuall y 
depends on the choices of the atom-photon coupling Q.^ and the scanning frequency step Awc, as seen in Eqs. (|^ 



(5.22) and (6.23). It is noticed in Eq. (8.4) that the signal rate is proportional to the number of axions contained 



in the conversion cavity {palTna)Vi. The dependence of the signal rate on the axion mass is further specified by 
considering the relation for the conversion cavity volume 



(8.5) 



This relation is due to the fact that the diameter of the cavity is proportion al to wh ile its lengt h is fixed. 

To summarize these arguments, the optimal situation, as given in Eqs. ( 6.21 ), (|6.22 ) and ( 3.23 ), can be set up 
experimentally for the dark matter axion search with Rydberg atom cavity detector. Then, the counting rates are 
expected to behave with respect to the changes of the relevant parameters as 



Rs 

Rn 



VCajjB^ffQPa ma Pa, 



i.7) 



Here, the noise rate is proportional to the thermal photon number nc[ma/Tc] which is determined as a function of 
the ratio nia/Tc. The resonant value a^a ~ (70/7)^^ ^ Pa^Q~^ is obtained for the axion-photon-atom interaction 
with the small enough axion detuning lAtJaj < Acjc/2 ~ 7a. (This is valid as long as the quantum uncertainty of 
energy ~ ?i/itr due to the finite atomic transit time is smaller than the axion width ?i7a, as ensured in the actual 
detection system.) It should also be remarked that the small detuning Aw;, of the atomic transition frequency only 
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shifts slightly the location of the peak of the axion signal from Aa'a(peak) = to Ati;a(peak) ~ Aw;,. By taking these 
signal and noise rates in the optimal case, the measurement time Ai at each scanning step is estimated with respect 
to the relevant parameters as 



B-^Q-^Plmlp-^ {R,/R^ > 1) 



The total scanning time is then estimated with the appropriate scanning step Acuc ~ 7a as 

ttot ~ At/Pl (8.9) 

Here the negative power of (3a appearing in the right-hand side is indeed compensated by the positive power of Pa 
contained in Eq. ( ^.8[ ) for At. The sensitivity seems to be improved apparently for the high atomic velocity in Eqs. 
( ^.8[ ) and ( |8.9| ). It should, however, be remarked that the condition h/tti- < hja for the resonant axion-photon-atom 
interaction is no longer ensured if the atomic velo city beccomes too high. The high atomic beam intensity is even 



required in such a case in accordance with Eq. (8^). Hence, we will in fact see in the next section that the preferable 
atomic velocity is w '--^ lOOms""'^ — lOOOms"^ for the actual experimental apparatus. 

Some essential features have been discussed so far concerning the detection sensitivity of the Rydberg atom cavity 
detector. In the next section, they will indeed be confirmed by detailed numerical calculations for the realistic situation 
with continuous atomic beam. 



IX. NUMERICAL ANALYSIS 



Numerical calculations have been performed to determine precisely the quantum evolution of the axion-photon-atom 
system in the resonant cavity, where some practical values are taken for the experimental parameters such as 

ma = 3x IQ-^eY - 3 x IQ-^eV, 

= 3 X 10"'' - 3 X 10-3, 
Pa = /Ohaio = 0.3GeVcm-3, 
Tc = lOmK - 15mK, 
Q = 1 X 10^ - 7 X 10"^, Seff 4T, 
Vi = SOOOcm^ {ma/lQ-^eVy^ , L = 0.2m, 

/Ryd = lO^S-l - lO^S-l, 

V = 350ms-i - lOOOOms-i, 
r2 = 5 X lO^s-i, n = IQ-^s. 



The steady state is realized by solving the master equation ( 5.18| ) for a long enough time interval MSt {M ^ 1). It 



in practice appears independently of the choice of initial value A/'(0) due to the finite damping rate s of t he axions, 
photons and atoms. The spatial distribution of the excited atoms in the cavity is calculated by Eq. ( 7.16| ) with this 
steady solution. The contributions of the galactic axions and the thermal photons p\j\x) are depicted together 

m Figs. I and | for the case of DFSZ axion with the tanh-type and sine-type electric field profiles, respectively, 
where the relevant parameters are taken as ttIq — 10-^eV (wc — ^a), Tc = 12mK, Q = 2 x lO'', = 0-1 ^'^d 

V = 350ms-^. These results obtained with K ^ \Q and M = 20 in fact appear to be smooth enough indicating that 
the steady state in the case with continuous atomic beam is realized in a good approximation. It is here remarkable 
that the form of electric field profile is reflected in the distribution of the excited atoms. In the following calculations, 
we take the tanh-type electric field profile for definiteness. The essential results concerning the interaction between 
axions and atoms mediated by photons in the resonant cavity, however, hold even in the cases with more realistic 
electric profiles. 

We can compute the signal and noise rates with the steady solutions obtained in this way. (The following calculations 
are made for if = 5 and M = 10 with the tanh-type electric field profile. The results are changed less than 10% even 
if larger K and M are taken.) It is readily checked by these calculations that the signal rate Rg is indeed proportional 
to c^^^i?gfjVi(pa/ma), as shown in Eq. (^^). Hence, the signal rate for the KSVZ axion is larger than that for the 
DFSZ axion by a factor c2^^(KSVZ)/c2^^(DFSZ) ~ 7. 

In Fig. ^, we plot the signal and noise rates, Rg and Rn, together depending on the atomic beam intensity /Ryd (and 
the corresponding atom-photon coupling fJjv), where the relevant parameters are taken as iria = lO^^eV (ojc = Wa), 
= 2 X lO'', L = 0.2m, v = 350ms-\ n = 5x lO^s'^ and Tc = lOmK, 12mK, 15mK. The noise rate is clearly 
proportional to the thermal photon number nc[ma/Tc], which is fic — 8.9 x 10"^, 6.2 x 10"'^, 4.3 x 10"^ at Tc = lOmK, 
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12mK, 15mK, respectively, for the axion mass iria = 10~^eV. In the weak beara intensity region, the signal and noise 
rates increase almost proportional to Inyd = N/ttr- On the other hand, for sufficiently strong beam intensities the 
noise rate is saturated to certain asymptotic value. In this case, the atoms interact reversibly with the photons in the 
cavity so that the equilibrium with nb[c — *■ &] ~ is realized. (In practice, due to the finite atomic lifetime the total 

number of excited atoms nb[c —^b]= p^f^\x)dx obtained from the thermal photons is slightly different from fie.) 
As for the signal rate , we clearly observe in Fig. ^that it is optimized for certain atomic beam intensity corresponding 



to the condition (8_^) for the atom-photon coupling. If the atomic beam intensity is too strong, the signal rate even 
decreases due to the Rabi splitting. It is hence important to adjust the beam intensity so as to attain the optimal 
condition for the signal. 

It is also checked that the lines representing the noise rate Rg versus beam intensity /Ryd shifts horizontally by 
changes of the intrinsic atom-photon coupling fl. Although the theoretical estimate of the electric dipole transition 



matrix element d is available ^0|] for calculating fl with Eq. (4.12), some ambiguities would be present in such a naive 
theoretical estimate of fl. Hence, it may rather be necessary to determine Q experimentally by fitting the thermal 
noise data in the weak beam intensity region to the expected lines, as shown in Fig. ^ for Q — 5 x lO'^s^^, which are 
calculated by varying ft around some plausible value suggested by the theoretical calculation of d . 

Another aspect should be pointed out concerning the counting rates of the excited atoms. As given in Eqs. d?!^ ) 



and ( 7.18 ), the noise rate is calculated by the formula i?„ = vp^^\L) with the excited atom density around the exit 
of cavity {x — L). Then, it may be noticed in Fig. || that the saturated value of i?„ in the strong beam intensity 
region, e.g., 0.23s^^ for T^. = 12mK, is in fact significantly larger than the naive estimate nc/itn e.g., O.lls"^ for 
nc[10~^eV/12mK] = 6.2 x 10^^, v = 350ms^^ and L — 0.2m. This enhancement of the counting rates, which is 

expected for the signal as well as the noise, is due to the fact that the densities of excited atoms pj,"' {L) and pj^' (L) 
around the exit of cavity become higher than the average values, as seen in Figs. ^ and ^. The nonuniform density 
of the excited atoms in the cavity is indeed brought by the fact that the atoms in the continuous beam arc passing 
through the spatially varying electric field with finite transit time. 

The behavior of the signal rate Rs with respect to the axion detunig Awa is shown in Fig. |^. Here some typical 
values are taken for the atom-photon coupling, ^n/i — 0.1, 0.3, 0.7, and the other relevant parameters are chosen 
as nria — 10~^eV, (3a = 10~^, Q ~ 2 x W^, L — 0.2m and v = 350ms^^. This behavior almost agrees with that of 
the form factor in Fig. ^ which is obtained for the simple case. (Note that the log-scale is taken for Rg in Fig. |^.) 
It is therefore confirmed by these numerical calculations that the resonant axion-photon-atom interaction takes place 
if the axion detuning is small enough as lAcJaj 7^ with the optimal atom-photon coupling ^n/i ^ {la/lY^'^ (as 
long as h/tf[ < h'ya)- It is also checked that if there is a small detuning Acuh of the atomic transition frequency, this 
resonant condition is slightly modified as \Au)a — AoJb\ — \uja ~ uJb\ ^ 7a- 

As clearly observed in Fig. ^, a salient feature which is found by these calculations for the case with continous 
atomic beam is the appearance of a ripple structure in R^ versus Awq with the relatively strong atom-photon coupling, 
^n/i > 0.1 for the present choice of parameters. It is realized that this structure is brought from the narrow axion 
width 7a P'^rria/h. We have checked that this fine structure indeed disappears if the axion spectrum spreads much 
wider with higher mean velocity such as /?a = 3 x lO^-^. 

The signal rate Rg versus axion detuning AcJq is calculated for some typical axion velocities, /3a = 3 x 10"'*, 1 x 
10-3,3 X 10-3. The results are shown in Fig. |, where ma = 10-^eV, Q = 2 x 10", 17^,/7 = 0.1, L = 0.2m and 
V = 350ms^^ are taken. It is here clearly observed that this signal form factor exhibits the structure determined 
by the energy spread of axions ?i7a ~ /3^ma. The galactic axion spectrum can actually have some narrow peaks, as 
pointed out in pl[ |. Then, by taking a small enough scanning frequency step Aojc, such peaks may be observed in the 
present detection scheme as well as the heterodyne method . The energy resolution is now limited by the qauntum 
uncertainty ~ h/ttr, which could be improved by utilizing the lower velocity atomic beam. 

The dependence of the axion signal Rs on the atomic velocity v is also shown in Fig. ||for nia — lO^^eV, f3a = 10^^, 
Hn/j — 0.1, L = 0.2m and v = 350ms-^, 2000ms-^, lOOOOms-^. Here, we note that for high atomic velocities such 
as V = lOOOOms-^ this axion signal form factor has a width much larger than that of axions. It is in fact determined 
by the energy uncertainty ~ h/tn (e.g., ~ 0.1?i7 for v ~ lOOOOms^^) due to the atomic motion. 

The measurement time At required for 3f7 level at each scanning step is shown in Fig. |^ depending on the atomic 
beam intensity /Ryd. Here the axion detuning is taken for definiteness as |Awa| < Aujc/2 with the scanning frequency 
step 



which is of the order of axion width 7a for /3a 10 ^ meeting the resonant condition ( |6.23 ) for the axion signal 



Then, the total scanning time ttot is estimated, as shown in Fig. Ml These calculations are made for the DFSZ axion 
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by taking some typical values of Tc, Q and v, which are indicated in Figs. ^ and|To[ The other relevant parameters 
are chosen as rUa = 10~^eV, pa = 10""^, ll = 0.1, L = 0.2m and v = 350ms^ . Here, we can see that t he b eam 
intensity to optimize the detection sensitivity changes with respect to Q and v according to the relation (8.3), as 



discussed before. The optimal estimates of the measurement time and the total scanning time also show roughly the 



dependence on Q and v as indicated in (3.8) and (8.9). For the KSVZ axion, the sensitivity is much better by a factor 

~ 72 ~ 50 {Rs/Rn < 1) or ~ 7 (Rs/Rn ^). 

The detection sensitivity for Bcr is finally estimated over the axion mass range nia = 3 x 10"^eV — 3 x 10~^eV, which 
can be searched with the present type of Rydberg atom cavity detector. In these calculations, the optimal atomic 
beam intensity is taken as Inyd = 4 x 10^s~^(TOa/10~^eV)2(Q/3 x 10^)^^ in accordance with the relation (pT^), which 



is also read from Figs. || and |l^. The volume of the conversion cavity is taken as Vi — 5000cm'^(r7i(i/10 ^V) ^ by 



considering the relation (8.5). The other relevant parameters are chosen as f3a — 10~'^, L = 0.2m, v = 350ms^^, 
Awc = 0.057 s-ncl Tc = lOmK, 12mK, 15mK. The results are shown in Figs. |ll| and [l^. It should here be noted that 
the Q factor can actually increase for lower frequencies. In order to take into account this property of the Q factor, 
we have assumed for example a relation 

Q{ma) = 3 X 10^(10-^eV/m„)2/3, (9.2) 

providing Q{ma) ^ 1.4 x 10^ - 6.7 x 10'' for 771^ = 3 x lO^'^gy ~ 3 x lO^^eV. The results obtained by taking a fixed 
value and this relation for the Q factor are plotted toghether by solid and dotted lines, respectively, in Figs. ^ and 



Here we clearly find that if the Q factor increases as given in Eq. (9.2), the detection sensitivity can be improved 
significantly for the lower axion masses nia ^ 10~^eV. It is also observed in most cases that the the measurement 
time At and the total scanning time itot once become local minimum for certain axion mass around nia = 10~^eV. 
As the axion mass gets smaller from this minimum point. At and itot increase until reaching a local maximum, and 
then turn to decrease again. This feature is unders tood by noting the behaviors of the signal and noise rates with 
respect to the axion mass. In fact, as shown in Eq. ( ^.6| ), the singal rate increases approximately as Rg ~ (with 
fixed Q) when the axion mass gets smaller. This increase of the signal rate is eventually overwhelmed by the more 
rapid increase of the noise rate i?„ ~ g-ma/feBTc Jqj. raa/k^Tc S> 1 proportional to the thermal photon number fie. 
Then, as the axion mass gets smaller to be comparable to the cavity temperature Tc, the increase of the noise rate 
becomes gradually moderate, and the singnal rate begins to dominate again. 

We can now conclude with these detailed numerical calculations. If the galactic axion search is made by utilizing 
the Rydberg atom cavity detecter, the DFSZ axion limit can be reached in frequency ranges of 10% around the axion 
mass TTLa ~ lO^^eV — lO^^eV for reasonable scanning times. The optimal condition for the detection sensitivity is 
attained by cooling the cavity down to a temperature Tc ~ lOmK and adjusting the experimental parameters such as 
the atomic beam intensity iKyd and velocity v and also the scanning frequency step Aujc- 



X. SUMMARY 



We have developed quantum theoretical calculations on the dynamical system consisting of the cosmic axions, 
photons and Rydberg atoms which are interacting in the resonant cavity. The time evolution is determined for 
the number of Rydberg atoms in the upper state which are excited by absorbing the axion-converted and thermal 
background photons. The calculations are made, in particular, by taking into account the actual experimental situation 
such as the motion and uniform distribution of the Rydberg atoms in the incident beam and also the spatial variation 
of the electric field in the cavity. Some essential aspects on the axion-photon-atom interaction in the resonant cavity 
are clarified by these detailed calculations. Then, by using these results the detection sensitivity of the Rydberg atom 
cavity detector is estimated properly. The present quantum analysis clearly shows that the Rydberg atom cavity 
detector is quite efficient for the dark matter axion search. 
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APPENDIX A: ANALYTIC SOLUTION 



Basic properties of the axion-photon-atom interaction in the cavity can be examined by considering the simple 
case where the atom-photon coupling does not depend on the time (i.e., the spatial variation of the electric field in 
the cavity is not considered). In this case, the analytic solution is obtained for the particle number matrix Afij{t) as 
follows 0. 

We first consider a more general case with K atomic bunches in the incident beam where the atom-photon couplings 
Qi - fix, which may be defferent each other, are independent of the time. The effective Hamiltonian Ti. given in Eq. 
(7.2) with constant f2i - fix is diagonalized by a nonsingular linear transformation T (not unitary for nonhermitian 

ny. 



W = T-iWT = diag.(Ai,A2,. 
Accordingly, the relevant matrices are introduced by 



(Al) 



(A2) 
(A3) 



where M'{t) = U-'^*{t)Af{t)ir_L^ (t), as given in Eq. ( |5.19D , with U{t) = exp[~int\ for the time-independent Ti. 
Then, the master equation ( |5.18| ) is reduced to a set of unmixed equations, 



dt 



These equations are readily solved with an appropriate initial value A/'(0) = N'{Q) as 
where 



A* - A, 



-i(A*-Aj)t 



1 



By using this solution for Af'{t), the original Af{t) is caluclated as 

U{t)=U*{t)[U{0)+C{t)]U^{t), 

where 

C{t) ^T*Cit)T^, 

U{t) ^ T exp[-im]T-\ 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 
(A9) 



The above procedure with Eqs. (A5) - (A9) to find the solution may be viewed as an inhomogeneous linear transfor- 
mation from Af{0) to JV{t) depending on the form of H: 



(AlO) 



The analytic solution obtained in this way can be used in a good approximation for a sufficiently small time interval 
where the variation of the atom-photon coupling is neglected. (This is equivalent to approximate the electric field 
profile f{x) wit h a re levant step-wise function, if the time interval is taken to be St = Sx/v.) Then, by applying the 
transformation (AlO) for a time sequence t = t^+i > ty > t^-i > ■ ■ ■ > Iq with small intervals — ifc-i (1 < k < z + 1) 
the approximate solution for M(t) may be obtained as 



jvito)"'^'^ mi) 

We now consider the specific case of = 1 with il{t) = iljv, where the Hamiltonian is given by 

ujb - ^7b 

H = I r^AT Uc - |7c K. 

K U!a- |7a 



H(ti 



u 



(All) 



(A12) 
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Then, the transformation matrix U {t) is calculated from Eq. ( [A9| ) as 

^..w-E4e"^''=*' (A13) 



with 



g1,=7l,T^\ (A14) 

where the sum is not taken over the index k. This transformation matrix may also be expressed in terms of the 
Laplace transform C |lj] as 

U{t) = C-^[{sl + in)-^]. (A15) 
The coefheients gfj are readily found by using the Laplace transform as 

4= \im{s + tXk){sl + tn)~j\ (A16) 

where in the leading orders of k 

SaSc -iflNSa -K^N 

A"^(s) ( -iflNSa SaSb -IKSh | (A17) 



with 



Si = s + iuJi + ^74, (A18) 



and 

A(s) = det(sl + in) 

= (s + iXi){s + iX2){s + iXs). (A19) 

The respective particle numbers are then obtained with the initial value M{0) — diag. (0, nc, na) as 

Uiit) = Mnit) = ricit)nc + na{t)na, (A20) 

where 



with 



1 - 



(A21) 



A™„ = -i(A;;-A„). (A22) 

The eigenmodes of TL are readily found as follows in the limit of k ^ 0, which is indeed sufficient for K/7 < 10~^^ 
or so: 

Xi^LOb- '-lb - Abe + [{Abcf + , (A23) 

X2=u;c- ^7c + Abe - [{Abc f + n%] , (A24) 

X3^UJa-^-/a, (A25) 

where Abe = ^{^b ~ ^c) ~ \(rib ~ 7c), and z^/^ = exp[z arg(z)/2]. The eigenmodes Ai and A2, which mainly 
consist of the atom and photon, are determined by the equation 



A - + -76 ) [X-ujc + 7;lc] -^'n^O- (A26) 



The eignemode A3 is, on the other hand, almost identical with the axion mode. These eigenmodes of Ti. satisfy the 
conditions 

Re[A,„„] = -(Im[A™] + Im[A„]) > 0, (A27) 
so that the respective particle numbers ni{t) converge to certain asymptotic values for t 00. 
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APPENDIX B: L ANGEVIN PICTURE 



We here describe the quantum evolution of tlic axion-photon-atom system in the Langevin picture p4[ |, which is 
supplementary to the treatment in the Liouville picture reproducing the same results for the quantum averages of the 
relevant quantities. The Langevin equation for the damped oscillators is given by 



(Bl) 



The external forces Fi are introduced for the Liouvillian relaxations. These external operators obey the following 
correlation property |ig| ]: 



(i^/(T)^;(r')) =%7.^.<5(r-T'), 
(^^;(r)g,(0))^0(r,r'>0). 



(B2) 



The solution of the Langevin equation ( p3l|) is readily found with the transformation matrix U (t) representing the 
time evolution due to Ti{t), as given in Eq. (^.20): 



q,it)=U,,it)q,iO) + / [U{t)l('\T)]^^F,iT)dT. 



(B3) 



Then, the time evolution of the particle numbers is calculated with this solution (B3) of the Langevin equation and 
the correlations (B2) as 



Af{0) 



(B4) 



kl 



where Dij = ^riiSij. This solution for JV{t) just coincides with that given in Eq. ( ^.24 ) which is obtained in the 
Liouville picture. 
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FIG. 1. A schematic diagram of the Rydberg atom cavity detector. 
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FIG. 2. The form factors for the axion conversion versus the cixion detuning for the case with constant atom-photon coupUng. 
The axion width is taken to be 7a = O.O27 with = 10~^ and Q = 2 x 10*. 
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FIG. 3. The distribution of excited atoms in the cavity for the tanh-type electric field profile f{x). The relevant parameters 
are taken as rua = 10~®eV {ojc = uja), Tc = 12mK, Q = 2 x 10*, f2jv/7 = 0.1 and v = 350ms~^. The contribution of axions 
(x) and that of thermal photons pj/' (a;) are shown together. 
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FIG. 4. The distribution of excited atoms in the cavity for the sine-type electric field profile f{x). The relevant parameters 
are taken as rua = 10~®eV {ojc = uja), Tc = 12mK, Q = 2 x 10*, f2jv/7 = 0.1 and v = 350ms~^. The contribution of axions 
(x) and that of thermal photons pj/' (a;) are shown together. 
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FIG. 5. The signal and noise rates versus the atomic beam intensity (atom-photon couphng). The relevant parameters axe 
taken as rua = 10"^eV (wc = Wa), Q = 2 x 10*, L = 0.2m, v = 350ms"^ and Q = 5x lO^s"^ 
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FIG. 7. The counting rate of signal versus the axion detuning for various values of the axion velocity Pa- The relevant 
parameters are taken as rua = 10~®eV, Q = 2 x 10*, Qn/j = 0.1, L = 0.2m and v = 350ms~^. 
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FIG. 8. The counting rate of signal versus the a^xion detuning for various values of the atomic velocity v. The relevant 
parameters are taken as ma = 10~®eV, j3a = 10~^, Q = 2 x 10*, Ojv/t = 0.1 and L = 0.2m. 
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DFSZ axion ma=10 eV 
v=350ms"^ Q=2xl0^ (solid) 
v=350ms"^ Q=4xl0^ (dotted) 
v=2000ms"^ Q=2xl0^ (dashed) 
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FIG. 9. The measurment time for 3a at each frequency step depending on the atomic beam intensity, which is estimated in 
the case of DFSZ axion. Here, some typical values are taken for Tc, Q and v, and the other relevant parameters are chosen as 
rua = IQ-^eV, /3a = 10"^, fiiv/l = 0.1, L = 0.2m and v = SSOms'^ 
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DFSZ axion ma=10 eV 
v=350ms"^ Q=2xl0^ (solid) 
v=350ms"^ Q=4xl0^ (dotted) 
v=2000ms"^ Q=2xl0^ (dashed) 
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FIG. 10. The total scanning time {3a level) for the DFSZ axion depending on the atomic beam intensity. Here, some typical 
values are taken for Tc, Q and v, and the other relevant parameters are chosen as rua = lO-^eV, /3„ = 10-^ njv/7 = 0.1, 
L = 0.2m and v = 350ms~^. 



33 



10' - 



< 



10" - 



10^ 



■ 1 

— 


1 1 1 1 1 i y 

n;^ / / — 

/ - 




DFSZ axion v=350ms~^- 




Q=3xlO^ (solid) 




Q(ma) (dotted) 




Tc=10mK (lower) 




12mK (middle) 


1 


15mK (upper) 
, , 1 1 



^-5 



10" 
ma(eV) 



FIG. 11. The measurment time for 3a at each frequency step depending on the cixion mass, which is estimated in the case 
of DFSZ axion. The Q factor is fixed (soUd lines) or varies with the axion mass (dotted lines). 
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FIG. 12. The total scanning time (3cr level) for the DFSZ axion depending on the axion mass. The Q factor is fixed (solid 
lines) or varies with the axion mass (dotted lines). 
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